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Gravitational waves imprint apparent Doppler shifts on the frequency of photons propagating 
between an emitter and detector of light. This forms the basis of a method to detect gravitational 
waves using Doppler velocimetry between pairs of satellites. Such detectors, operating in the milli- 
hertz gravitational frequency band, could lead to the direct detection of gravitational waves. The 
crucial component in such a detector is the frequency standard on board the emitting and receiving 
satellites. We point out that recent developments in atomic frequency standards have led to devices 
that are approaching the sensitivity required to detect gravitational waves from astrophysically 
interesting sources. The sensitivity of satellites equipped with optical frequency standards for Doppler 
velocimetry is examined, and a design for a robust, space-capable optical frequency standard is 
presented. 


I. INTRODUCTION 

The detection of gravitational waves (GWs) is one of the prime goals of gravitational physics today. Although there 
is tantalizing indirect evidence for the existence of gravitational waves [1], they have never been detected directly. At 
present, there are two working approaches to detecting GWs - pulsar timing arrays (PTA) in the nanohertz GW 
frequency band [5], and terrestrial optical interferometers in the dekahertz band [3]. These approaches, however, are 
insensitive to the astrophysically natural /rHz-mHz band of GW frequencies. PTAs are insensitive in this band due to 
the the long transit time between the pulsar and earth, and the need to average pulsar signals for a long duration to 
obtain adequate signal-to-noise ratios. Terrestrial interferometers, on the other hand, are affected by low-frequency 
noise sources which limit their operation below ^ 1 Hz. Upcoming missions such as eLISA [3], building upon the LISA 
concept [315] of measuring the distance between distant inertial test masses using ultra-stable lasers, are expected to 
probe the mHz GW band. Other proposed approaches to detecting GWs in this frequency band include the use of 
very large atom interferometers |5]- 

A complementary approach to probing the /rHz GW band is to use satellite Doppler velocimetry (SDV). This 
approach derives from the pioneering analyses of Kaufmann |7] , and Estabrook and Wahlquist |8] . We refer the reader 
to Smarr et al. [3 for a description of one possible detection concept, and nni for a comprehensive recent review. 
Kaufmann’s original analysis |7| showed that a GW, interacting with light propagating from an emitter to a detector, 
imprints a small apparent Doppler shift on the light wave. This apparent Doppler shift is proportional to the difference 
in GW amplitudes at the emitter and detector (times an order unity angular factor). If this small frequency shift can 
be resolved, it offers an attractive approach to detecting GWs. The apparent Doppler shift imprinted by the GW is 
oscillatory at the GW frequency, and therefore it is necessary to resolve the small Doppler shift on a timescale less 
than the GW period so that the signal is not averaged away to zero. This places strong constraints on the performance 
of the frequency standards carried by the emitter and detector, requiring them to reach a high level of stability within 
a relatively short averaging time (1000s of seconds). 

The purpose of this article is to point out that the demonstrated performance of present-day optical frequency 
standards makes the SDV approach competitive with complementary detection proposals for low-frequency gravitational 
waves. In the following, we first briefly review the physics of SDV, followed by an analysis of the performance achievable 
with optical frequency standards. We then present a design for a portable optical frequency standard that can be used 
as a component of SDV detectors of gravitational waves. 


II. GRAVITATIONAL WAVE DETECTION USING SATELLITE DOPPLER VELOCIMETRY 

Consider a gravitational wave with a wave-vector k propagating along the 0 -axis in a region of space. The metric in 
this region is 


= dt^ — {dz^ + dx^ + dy"^) — h+[fc (0 — t)]{dx^ — dy^) — 2hx [k{z — t)]dxdy, (1) 

written in the transverse-traceless gauge. [35] Here and /ix denote the amplitudes of the two polarizations of the 
GW. The cyclic coordinates of the action integral of a free particle (or photon) in this region are x,y,z + t, leading to 
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the conservation of the momentum components px.,Py,Pz +Pt- In addition, a photon’s four-momentum satisfies the 
null constraint p ■ p = 0.[27] 

For a photon propagating with constant components Px = a, Py = I3,pz + Pt =7, the null condition p-p = 0 implies 
that its four-momentum has the form 


p = 




( 2 ) 


where 


_ a'^ + + (a^ — /3'^)h^ -|- 2a/3/ix 

^2 

= tan^ 6* (1 -I- cos2(j)h+ + sin2(j)hx) 


(3) 


with the constants a = po sind cos (/),/? = po sin d sin 0,7 = pgcos^ written in terms of spherical polar parameters 

po,6>,0. 

The photon’s frequency measured by a free-falling emitter (or receiver) v = i^p ■ (1; 0, 0, 0) = ^(1 -I- S). The ratio 
of the frequencies emitted by the emitter e, and detected at the receiver r, is therefore 

— = « 1 -|- sin^ 0 (cos 20 -|- sin 20 Ahx ), (4) 

Vr 1 + Or 

where Ah^ = /i+(e) — /i+(r), Ahx = (e) — hx (r) are the differences between the GW strain amplitudes at the emitter 

and receiver. The frequency ratio r'e/r'r is different from unity when or Ahx is non-zero. The GW amplitude 
appears as an apparent Doppler shift at the receiver. It can be measured by precisely comparing the frequency of light 
received from the emitter to that of a stable frequency standard at the receiver. Therefore the fundamental limit to 
the performance of such a GW detector is the stability of the frequency standard carried by the emitter and receiver, 
evaluated on a timescale corresponding to the period of the GW. Figure shows the single-cycle GW sensitivity of 
an emitter-receiver pair on board satellites in solar orbit (here onwards a “GW detector”), for various values of the 
stability of the frequency standard. The performance of a SDV detector has a slower low-frequency roll-off compared 
to laser phase-measurement-based detection methods, and therefore offers a complementary method for low-frequency 
GW physics. (For the parameters chosen in Figure both the SDV response curves are more sensitive than the 
proposed eLISA detector [i] in the 1-100 /rHz GW band.) 

As was pointed out by Kaufmann [7], and Estabrook and Wahlquist [ 8 ], note that a detectable frequency shift is 
obtained in spite of the fact that there is no GW perturbation on the local clock rates at the emitter and receiver 
(the t-component of the metric is unaffected by the GW). The imprinting of the GW’s amplitude (which appears in 
the spatial part of the metric) onto the frequency of the photon is a consequence of the propagation of the massless 
photon through the perturbed spacetime. 

We note that the stability of the frequency standard and its high signal-to-noise ratio (SNR) are crucial, not its 
absolute accuracy. Further, well-characterized long-term drifts may also be tolerated, as the GW signal appears as an 
oscillatory apparent Doppler shift at a characteristic Fourier frequency within the sensitivity band of the GW detector. 
This opens up the possibility of using frequency standards that may not necessarily be the best for time-keeping 
applications, but instead trade absolute accuracy for portability, stability and SNR. 

The satellites will experience significant gravitational redshifts from being at different locations in the solar system’s 
gravitational potential (making their absolute accuracy moot anyway), which will lead to slowly varying frequency 
shifts. In addition, effects such as perturbations from planets will lead to time-varying frequency shifts due to their 
effect on the satellite trajectories. However, these effects will typically appear on slower timescales (~months for 
orbital perturbations) compared to the ^hour oscillation period of GW in the /rHz-mHz band. These spurious effects 
can therefore be clearly distinguished from GW signals. 

The primary source of noise in the Doppler signal arises from uncontrolled motion of the satellite platform (with an 
approximately white velocity spectral density of ^ 10“® m/s/v^Hz at /rHz-mHz frequencies), due to the effects of solar 
wind, radiation pressure, and the spacecraft’s mechanical imperfections m- We suggest a method to compensate for 
this effect, using the fact that the resultant displacement noise is enhanced by a factor oc 1 //: the displacement of an 
isolated free-falling test mass on board the satellite can be monitored, and the displacement record communicated 
between the satellites. This record is processed to correct the measured Doppler signals against the fluctuations in the 
relative velocity of the satellite platforms. By monitoring the displacement of the on-board inertial mass e.g. capacitive 
or optical interferometric sensing, the velocity noise can be compensated. With a displacement measurement sensitivity 
of 0.5 pm/Vilz (2 orders of magnitude less critical than the demonstrated performance of the drag-free control system 
in the LISA Pathfinder, in a similar frequency band m), the velocity noise can be compensated to < 10 ® m/s/-\/Hz 
at 100 /iHz. At this level, the residual platform velocity noise is no longer a significant concern for SDV detection, as 
shown in Figure 
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FIG. 1: Calculated response of a SDV GW detector, consisting of an emitter-receiver pair separated by 2 astronomical units 
(AU). The green solid line represents the single-cycle GW strain sensitivity (“burst sensitivity”) achievable with state-of-the-art 
atomic frequency standards m, and the blue dashed line represents the performance achievable with a frequency standard 
that has ten-fold improved stability - both lines are labeled by their SNR. Also shown is the limiting sensitivity due to residual 
velocity noise of the satellites, after corrections obtained from displacement sensing of an isolated inertial mass at the level of 0.5 
pm/y/Wz. The red line is the theoretically predicted characteristic strain for GW radiation from one class of potential sources 
showing that such sources are already detectable with state-of-the-art frequency standards. 


III. OPTICAL FREQUENCY STANDARDS 

Buoyed by the invention of the optical frequency comb m, optical frequency standards have made enormous leaps 
in performance over the last decade, and are presently the best frequency standards available [H] [15] . These devices 
generate a stable optical frequency by slaving a laser oscillator to a high-quality-factor atomic resonance. The current 
state-of-the-art is achieved by optical lattice clocks, where the reference atoms are cooled and trapped in standing 
waves of light El- Optical frequency standards are tremendously attractive as components of a SDV GW detector. 
Figure [^illustrates one possible method for performing frequency comparisons between the emitter and receiver in a 
GW detector equipped with optical frequency standards. 

At the emitter, a frequency comb is locked to an optical frequency standard (laser oscillator stabilized to an atomic 
frequency reference). The microwave frequency, obtained by dividing down the optical frequency using the frequency 
comb [16], is transmitted over the inter-satellite link to the detector. At the receiver, the received microwaves are 
amplihed by a low-noise amplifier, multiplied up to the optical domain using a frequency comb, and the resulting 
optical frequency compared to the receiver’s frequency standard. The microwave link can be operated with multiple 
frequencies to cancel the effect of plasma dispersion unKiT]. Alternatively, it is foreseeable that a high signal-to-noise 
ratio optical link can be operated between satellites in the near future [HI HI] , allowing direct comparisons between 
optical signals from the emitter and receiver. 

Most of the components involved in this scheme are compatible with operation on board a satellite, using currently 
available technology. This includes miniaturized optical frequency combs that have been developed in recent years |20j . 
The notable exception is the frequency standard, as optical frequency standards capable of being used on satellites 
have yet to be developed. In the next section, we address the feasibility of building and operating such a frequency 
standard. 
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FIG. 2: An example block diagram of a satellite Doppler velocimetry GW detector, with a simple one-way single-frequency 
link. In the emitter, a microwave frequency is obtained from the optical frequency standard, and transmitted over the link to 
the detector. At the receiver, the received microwaves are amplified, multiplied up to the optical domain and compared against 
the receiver’s optical frequency standard. Multi-path links and multi-satellite links can be used to enhance the GW detection 
sensitivity m, and multiple-frequency links can be used to suppress the effects of plasma dispersion }10^ . 


A. Towards a satellite-compatible optical frequency standard 

Due to their size and complexity, the state-of-the-art architecture for optical frequency standards (optical lattice 
clocks) may not be the best option for satellite-borne frequency standards (although cf. the efforts to miniaturize 
optical lattice clocks for measurements on the International Space Station m)- Other more portable designs, such as 
single-ion clocks [22] or solid-state clocks (23], might be better suited to being operated on a satellite. One of the main 
sources of the complexity of state-of-the-art optical frequency standards arises due to the long cycle times involved in 
interrogating the atomic reference using the laser oscillator. A high-performance laser oscillator, often with a linewidth 
of 1 Hz, must be used to interrogate the atomic transition for many seconds and the atoms must be held in place for 
the duration of the interrogation. By relaxing the constraints on the atomic frequency reference and reducing the 
measurement cycle time, the engineering of the frequency standard can be simplified. In this section, we present an 
alternative design for a portable optical frequency standard that could be useful as a component in a SDV GW detector. 
Satellites carrying high-peformance frequency standards would also be useful in precisely measuring gravitational 
properties of the solar system (such as gravitational redshifts) and constraining models of modified gravity. 

We have studied the design of an atomic frequency standard using the 916 nm ^5'o — ^D 2 two-photon transition 
in neutral calcium atoms. This transition has a linewidth T ~80 Hz, which relaxes some of the constraints on the 
laser oscillator required to interrogate it, and leads to a faster cycle time. The relevant atomic states and transition 
wavelengths in calcium are shown in Figure [^ The choice of this atomic reference also means that the frequency 
standard can be built using a single-stage magneto-optical trap (MOT), without incurring the complexity of the 
multiple stages of laser cooling and trapping required in optical lattice clocks. In contrast to schemes using single 
photon excitation, the thermal motion of the atoms in the MOT does not lead to significant frequency shifts in this 
design, due to the Doppler-free and recoil-free nature of degenerate two-photon excitation. The measurement scheme 
is compatible with low-power operation using robust and inexpensive diode lasers and a simple vacuum apparatus, and 
could lead to a portable and maintenance-free atomic frequency standard. We estimate that a frequency stability 
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FIG. 3: The two-photon ^So D 2 clock transition in neutral calcium (916 nm, 327 THz) is an attractive candidate for operating 
a simple optical frequency standard. The relevant transitions are shown in (a), and all the required laser wavelengths can be 
obtained from compact laser diodes. The measurement sequence involves the following steps: (b) atoms are cooled and trapped 
in a magneto-optical trap (MOT) operating on the 423 nm ^ Pi transition, (c) they are excited out of the ^So state by 

two counter-propagating photons obtained from the 916 nm laser oscillator, and (d) the excited state fraction is detected using 
laser-induced fluorescence on the ^D 2 F 3 cycling transition at 488 nm. 


5vjv « 10 !i/t can be obtained with an integration time of r seconds, with commensurate control over systematic 

frequency shifts. 


IV. SUMMARY 

The high performance of modern optical frequency standards can be leveraged in order to detect gravitational 
waves using satellite Doppler velocimetry. The performance of such GW detectors is comparable to, and complements, 
proposed space-based interferometric detectors, but relies on a simpler approach based on demonstrated technologies. 
The feasibility of operating such detectors hinges on the availability of portable and maintenance-free optical frequency 
standards, and we have described one path to constructing such a device. The development of satellite Doppler 
velocimetry for GW detection could offer a unique window into the gravitational wave spectrum, and address important 
open questions in gravitational physics. 
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Note added: Recently, somewhat similar proposals have been reported [Ml I25j . 
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